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Abstract 

All fundamental physical principles are confirmed by numerous exper- 
iments and practical certainty and have unambiguous interpretation. But 
physics of stars is based on few measured effects only. It gives some free- 
dom for figments of the imagination. The goal of this paper is to compare 
two alternative astrophysical models with measurement data. 



1 Introduction 

All basic problems of physics are thoroughly examined. The total accuracy of 
their solutions is confirmed by numerous direct experiments and their practical 
using. Any alterative approach to these basic problems seems not to be rea- 
sonable. At the first sight it must be applicable to astrophysics too, because 
it is a part of physics. But the situation here ia quite different. It is accepted 
in astrophysical community to consider the basis of modern astrophysics as ab- 
solutely reliable and steady. But this basis was developed in the past, when 
there were no possibilities to check them by measuring stellar parameters and 
many branches of physics, like plasma physics, did not exist. It forces to revisit 
the basis of stellar physics intently. Nowadays there is possibility to check ba- 
sic astrophysical problems by means of comparison of theoretical models with 
measurement data. The astronomers are able to observe and to measure a few 
parameters of stellar radiation and stellar moving, which obviously depend on 
the state of stellar interior. In the first place there are parameters of following 
phenomena: 

(1) the apsidal rotation in binary stars 

(2) the spectral dependence of solar seismical oscillations. 

From this point of view one can consider also the measurement of solar 
neutrino flux as one of similar phenomenon. But its result can be interpreted 
ambiguously because there are a bad studied mechanism of their mutual con- 
versation and it seems prematurely to use this measurement for a stellar models 
checking. 
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2 Two models of stellar interior 



It is generally accepted to think that the equilibrium of substance inside a star 
can be described by the Euler equation: 

VP + 7 g = 0, (1) 

where g is the gravity acceleration, P is the pressure and 7 is the substance 
density. Here the density satisfies 

- iTrG-f = div g, (2) 

where G is the gravity constant. Based on this equation astrophysicists supposed 
that the pressure and the temperature inside a star are growing monotonic with 
the increase of the depth. As plasma inside a star can be considered as an ideal 
gas at the pressure P = n kT (n and T are its density and temperature), it 
is supposed by all astrophysical models that the temperature and the plasma 
density grow in direction to the star center. Here the temperature reaches tens 
millions Kelvin and the density is approximately hundred times greater than 
its averaged value. It is a fundamental statement of modern astrophysics. It is 
important to underline that although the matter inside a star can be described 
as an ideal gas, it is not a gas. It is electron-nuclear plasma. Historically 
the Euler equation (Eq.Q) was formulate and applied to the astrophysical 
objects at the time when the term "plasma" did not exist and basic concepts of 
plasma physics were not developed. The features of electron-nuclear plasma can 
not be described by this equation. The equilibrium state in general case must 
take into account the balance of all forces applied to the system. As particles 
of plasma have masses and electric charges, the gravity action on plasma can 
induce its electric polarization V . Taking into account the gravity induced 
electric polarization (GIEP), the equilibrium equation obtains the form: 

VP + 7g + 47T? divT = (3) 

Analyzing this equation, it is easy to see that there are at last two possibilities. 
At first, the plasma body can exist in self gravitating field in a state (which can 
be called Eulerian), when J 1 = and an equilibrium is determined by Eq.Q). 
But another equilibrium state is possible when the gravity force is balanced by 
the electric force: 

7 g + 4ttT divV = (4) 
and the pressure inside the plasma body is constant: 

VP = 0. (5) 

(The density of plasma under this condition must be constant too 7 = const). 
Since the polarisation which is non-uniformly distributed in space, can be con- 
sidered as a distribution of "bond" charges 

divf = Pbond, (6) 
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we can rewrite the equilibrium condition in the following way: 



7S + Pbond E = (7) 

where the intensity of the electric field can be expressed through the density of 
"bond" charge: 

^Pbond = div E. (8) 

and we can describe further the equilibrium state by Eq. J7J which looks as more 
convenient at consideration. 

The theoretical consideration based on the Thomas-Fermi approximation 2 
states that the hot dense plasma in a gravity field must have the equilibrium at 
the following value of the bond charge electric density: 

Pbond = G 1 ' 2 1 (9) 

and at the electric field intensity 

E = gT> (io) 



3 The equilibrium of plasma at GIEP effect 

3.1 The equilibrium density of hot electron- nuclear plasma 

The self-gravity of a celestial body tends to compress the plasma in its central 
region. In this compressing a plasma density reaches up to 10 26 particles / cm 3 . 
Such a density corresponds to a minimum of the plasma energy 0. At tem- 
peratures about tens million Kelvin, plasma can be considered as an ideal gas. 
In accordance with the ideal gas definition, plasma particles interactions are 
neglected in this approximation. To take into account these interactions in next 
step of approximation, one can see that two mechanisms of interaction play the 
most important role. 

1) Electrons are Fermi-particles, and they obey the Fermi-statistics. They 
can not occupy levels in the energetic distribution which another electrons have 
taken. The value of correction for this interaction is known, it is given in Landau- 
Lifshitz course pQ . This correction is proportional to the density of particles in 
the first power and it is positive, therefore when one takes it into account, the 
plasma energy in this approximation is greater than the ideal gas energy at the 
same density and temperature. 

2) Electrons and nuclei have uniform space distribution at very high temper- 
atures (when interactions can be neglected completely). The same correlation in 
the space distribution between plasma particles appears at finite temperatures. 
It can be described by introducing a correlation correction which takes into 
account electrostatic interaction between nuclei and electrons. This correction 
is considered in Landau-Lifshitz course yQ. This correction is proportional to 
particle density in 1/2 and it is negative, therefore taking it into account leads 
to the decreasing of the plasma energy in comparison with the ideal gas energy 
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at the same density and temperature. If both these corrections are taken into 
account, we can see that there is a minimum of plasma energy at its density [2]: 



16 (Z + l) 3 



(11) 



V = 



9n r! 



Where Z is the averaged charge number of nuclei from which plasma consists, 
r B = e 2 is Bohr radius. The existing of equilibrium plasma density rj is not 
caused by GIEP effect directly. The conclusion about equilibrium density of 
hot plasma can be deduced from the standard statistical theory and has gen- 
eral meaning. Taking into account GIEP effect one concludes that the plasma 
density is constant. For equilibrium state this density is obviously equal to 77 
because it corresponds to the energy minimum. 

The equilibrium state with the density 77 is inherent to the dense plasma at 
high temperature (about 10 7 K). This temperature is characteristic for central 
region of a star only. It leads to the conclusion that the core with density 77 is 
placed in central region of a star and outside of the core there is a region where 
VP =/= and the density and the temperature are change to values which are 
characteristic for the star surface. 

3.2 Another equilibrium parameters of star cores 

The constancy of the pressure (VP = 0), which is characteristic for star cores, 
needs the constancy of the temperature (VP = 0). The value of the equilibrium 
temperature can be extracted from the temperature dependence of the energy. 
The temperature dependence of the star has two branches with different slopes. 
At high temperature, when the energy of radiation has a important role, a star 
as a whole has a positive heat capacity: its energy increases with increasing of 
its temperature. At a relatively low temperature the heat capacity of a star is 
negative. Here the star energy decreases while its temperature increases. This 
behavior of stellar heat capacity is a well known fact, it is discussed in Landau- 
Lifshitz course pQ. Of course, the own heat capacity of the star substance in each 
small volume is positive. One obtains the negative heat capacity of a star as a 
whole at taking into account the gravitational interaction. There is a minimum 
of the energy placed between these two branches of the temperature dependence 
of the energy with positive and negative slopes. It determines the equilibrium 
temperature of the star core [2] : 



A similar argumentation gives a possibility to determine the equilibrium mass 
and the equilibrium radius of star core [5] ■ 



Tw2-10 7 (Z + l) K. 



(12) 
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4 The apsidal rotation in binary stars 



4.1 The rotation of close double stars at Eulerian distri- 
bution of substance 

The apsidal rotation (or periastron rotation) of close binary stars is a result 
of their non-Keplerian moving which originates from the non-spherical form of 
stars. This non-sphericity has been produced by a rotation of stars about their 
axes or by their mutual tidal effect. The second effect is usually smaller and 
can be neglected. The first and basic theory of this effect was developed by 
A.Cleirault in the beginning of XVIII century. Now this effect was measured 
for approximately 50 double stars. According to Clairault's theory the velocity 
of periastron rotation must be approximately in 100 times faster if a matter 
is uniformly distributed inside a star. Reversely, it would be absent if all star 
mass is concentrated in a star centrum. To reach an agreement between the 
measurement data and calculations, it is necessary to assume that the density 
of the substance grows in direction to the centrum of a star and here it runs 
up a value which is a hundred times greater than a mean density of the star. 
Just the same mass concentration of the stellar substance is supposed by all 
standard theories of a star interior. It has been usually considered as a proof of 
astrophysical models. But it can be considered as a qualitative argument. To 
obtain a qualitative agreement between theory and measurements, it is necessary 
to fit parameters of the stellar substance distribution in each case separately. 

4.2 The apsidal motion of close binary stars at taking into 
account the GIEP effect 

In the absence of rotation a star would have a spherical core. The rotation 
transforms it in a oblate ellipsoid. Its oblateness can be calculated [3] and the 
velocity of periastron rotation can be obtain according to Clairault formulas: 



where V is the period of the ellipsoidal rotation of stars, U is the period of 
the periastron rotation. The parameter T is the period depending on world 
constants only: 



and the parameter 



depends on chemical composition of star cores. There Zi and Ai are the charge 
and the mass number of nuclei which are composing the plasma of i-star. Hence 
the velocity of periastron rotation depends on the chemical composition of a star 




(13) 




(14) 
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Figure 1: The distribution of binary stars on the value of ("P /U)^ /T) 2 . 



only and it decreases rather sharply with the increasing of Z, so the periastron 
rotation of a pair consisting of heavy nuclei must be indistinguishable. The 
calculations of £ for a some light nuclei are shown in Tab.l. 



Table 1. 



starl 


star2 


6+6 


composed of 


composed of 




H 


H 


.25 


H 


D 


0.1875 


H 


He 


0.143 


H 


hn 


0.125 


D 


D 


0.125 


D 


He 


0.0815 


D 


hn 


0.0625 


He 


He 


0.037 


He 


hn 


0.0185 



Here the notation " hn" indicates that the second component of the couple 
consists of heavy elements or it is a dwarf. The distribution of close binary stars 
on the value of (V /M)(V /T) 2 is shown on Fig^in the logarithmic scale. 

All these data and references were given to us by Dr.K.F.Khaliullin (Stern- 
berg Astronomical Institute) and are cited in [3]. On Fig^the lines mark the 
values of parameters fEn.qi3[)) for different pairs of binary stars. It can be seen 
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that the calculated values of the periastron rotation for stars composed by light 
elements which is summarized in Table 1 are in a good agreement with separate 
peaks of the measurement data distribution. It is important to underline that 
these results were obtained without using of any fitting parameters (and they 
are accurate in this regard). Reversely, the conventional approach to this prob- 
lem based on Euler equation can not give any explanation for this measured 
distribution. 

5 The solar seismical oscillations 

5.1 The seismical oscillations at the Eulerian distribution 
of a substance 

The measurements show that the Sun surface is subjected to a seismic vi- 
bration. The most intensive oscillations have the period about five minutes and 
the wave length about 10 4 km or about hundredth part of the Sun radius. Their 
spectrum obtained by BISON collaboration is shown on Fig[2 

These oscillations are a superposition of a big number of different modes 
of resonant acoustic vibrations. It is supposed that acoustic waves propagate 
on different trajectories in the interior of the Sun and have multiple reflection 
from surface. With these reflections trajectories of same waves can be closed 
and as a result standing waves are forming. Specific features of spherical body 
oscillations are described by the expansion in series on spherical functions. These 
oscillations can have a different number of wave lengths on the radius of a sphere 
(n) and a different number of wave lengths on its surface which is determined 
by the l-th spherical harmonic. It is accepted to describe the sunny surface 
oscillation spectrum as the expansion in series 

v rdm ~ Av (n + l - + e ) - 1(1 + 1)D + mAv rot . (16) 

Where the last item is describing the effect of the Sun rotation and is small. 
The main contribution is given by the first item which creates a large splitting 
in the spectrum (FigJSJ 

Av = Un+1,1 - Vn.l- (17) 

The small splitting of spectrum (Fig|2J depends on the difference 

5 VI = v n ,i - Vn-\,l+2 ~ (4/ + 6)D . (18) 

A satisfactory agreement of these estimations and measurement data can be 
obtained at [Sj 

Av = 120 fj,Hz, e = 1.2, D = 1.5 [iHz, Av rot = lfiHz. (19) 

To obtain these values of parameters Az^o, to Dq from theoretical models is 
not possible. There are a lot of qualitative and quantitative assumptions used at 
a model construction and a direct calculation of spectral frequencies transforms 
into a unresolved complicated problem. 
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Figure 2: (a) The power spectrum of solar oscillation obtained by means of 
Doppler velocity measurement in light integrated over the solar disk. The data 
were obtained from the BISON network |3j . (6) An expanded view of a part of 
the frequency range. 
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5.2 The oscillation of the solar core at the GIEP effect 



Fig. 2 shows a central part of the whole spectrum of solar oscillations which was 
obtained at a very high frequency resolution. The whole spectrum of solar os- 
cillations was obtained at a little worse resolution in the frame of the programm 
SOHO/GOLF 0, and it is shown in FigEJ 

The existence of this spectrum forces to change the view at all problems 
of solar oscillations. The theoretical explanation of this spectrum must give 
answers as minimum on four questions : 

1. Why does the whole spectrum consist from a large number of equidistant 
spectral lines? 

2. Why does the central frequency of this spectrum T is equal approximately 
to « 3.23 mHz? 

3. Why does this spectrum splitting / is equal approximately to 67.5 fj,Hz? 

4. Why does the intensity of spectral lines decrease from the central line to 
the periphery? 

The answers to these questions can be obtained from the Sun core model 
based on the GIEP effect. According to this model, the Sun core has the high 
constant density (Eq.JTTJ) with radius R, mass M at temperature T, which all 
depend on Z and A/Z only 

5.2.1 The sound oscillation of the Sun core 

Since the solar core is compressible, the main mode of its vibration should be 
elastic sound oscillations of its radius with the conserved spherical form of the 
core at the frequency fl s « V SO und/^~ The detail calculation shows that the 
frequency of this mode \7\ 

is depending on a chemical composition of the core only. The same separate 
frequencies of this mode of the sound radial oscillation (T = H s /2tt) for cores 
with the different A/Z a.t Z = 1 Z = 2 are shown in the third column of Table 
2. 
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Figure 3: (a) The spectrum of solar oscillations obtained by SOHO/GOLF 
group [S] . (6) The theoretically obtained spectrum with using Eq. (|25[) at Z = 2 
A/Z = 5. 
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4 


2.85 
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4.5 


3.02 






2 


5 


3.19 


The Sun 


3.23 



Table 2. 



The measured frequencies of the surface oscillations of separate stars |S] are 
shown in the right side of Table 2. 

Comparing the calculated and measured frequencies, one can conclude that 
the Sun core must be composed in general by hcllium-10. This conclusion 
doesn't look so confusing because the pressure acting in core is running to 
10 18 dyn/cm 2 and it can induce the neutronization process in plasma which 
makes neutron-excess nuclei stable. At this chemical composition we have 

=3.19 mHz. (21) 

The good agreement with the measurement data gives a possibility to argue 
that the central frequency of solar oscillation is related to the radial oscillations 
of its core. 



5.2.2 "Phonon-like"low frequency oscillations 

Another mode of the solar core oscillations is related to the existence of the 
equilibrium core density r\ fEci lllfl . Any deviations of the density from this 
equilibrium value, for example which are caused by radial core oscillations, in- 
duce a mechanism of density oscillations around this equilibrium value with 
frequency 

f V¥ 2 4 3/2 \Gm p ] /A\ T I 4 ' 5 ) 1 / 2 

w„ = < -a 3/2 —^t 2 - ( — > Z + 1 \ , (22) 

" \VT0(3/2)7 I r% \ \z/[ J J V 1 

It gives at Z = 2 and A/Z = 5: 

fn = T L = 66.0 fiHz. (23) 
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These oscillations of plasma are like plionons in solid bodies. The excitation 
of this mechanism can induce oscillations on multiple frequencies kuj^. Their 
intensity must be k times weaker because a population of according levels in 
energetic distribution is reversely proportional to their energy khw^. As a result 
these low frequency oscillations form spectrum 



^ - sin(ku! n t) (24) 



k 

k=l 

If these low frequency oscillations f v are induced by the sound radial oscil- 
lation with frequency T ', they will modulate them. The radial displacements of 
the solar core surface form the spectrum 

uji ~ sin Vl s t ■ — sin kuj^t w ^ sin Vl s t + — sin (Sl s ± kcj^t, (25) 

k=0 k=l 

This calculated spectrum is shown on Fig[3h>. It is can been seen from this 
figure that the model of the Sun based on GIEP effect gives a possibility to 
explain all basic details of the measured spectrum of oscillations and to obtain 
answers on all four questions which was formulated above. It is important to 
underline that the quantitative agreement between the calculated spectrum and 
the measurement data was obtained without using of any fitting parameters and 
only taking into account its chemical composition (Z = 2 A/Z = 5). 

6 Another measurements which results are de- 
pending on a star interior construction 

6.1 Stellar masses 

The models of star interior based on the Euler equation does not give any 
possibilities to estimate values of star masses. The theory based on GIEP effect 
obtains a direct way for the star mass calculation 2 

"^■ 5e (S) 1,2 (^) 3,2 (l> 2 "^ 6 - 47Afc " (f> 2 - (26) 



Where Men = ( ) m p = 3.71 • 10 33 <? is the Chandrasekhar mass. 



3/2 

The frequencies of natural oscillations of the Sun show that the star core 
mass is approximately equal to 1/2 of the whole stellar mass. It gives us a way 
for the estimation of stellar masses and for the comparison with measurement 
results. The calculated values of star mass depend on one coefficient A/Z only 

El 

There is no way to determinate the core chemical composition of far stars, 
but some predications for it are possible. At first, there must be no stars which 
masses exceed the Sun mass more than one and a half orders, because it accords 
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Figure 4: The mass distribution of binary stars |5]. On abscissa, the logarithm 
of the star mass over the Sun mass is shown. Solid lines mark masses agree with 
selected values of A/Z from Ea. (|26|) 



to limiting mass for stars consisting from hydrogen with A/Z = 1. Secondly, 
though the neutronization process makes neutron-excess nuclei stable, there is 
no reason to suppose that stars with A/Z > 10 (and with mass in hundred times 
less than hydrogen stars) can exist. Thus, GIEP-theory predicts that the whole 
mass spectrum must be placed in the interval from 0.25 up to approximately 
25 solar masses. These predications are verified by measurements quite exactly. 
The mass distribution of binary stars is shown on FigQ] 8 . (Using of this data 
is caused by the fact that only the measurement of parameters of the binary 
star rotation gives the possibility to determine their masses with satisfactory 
accuracy). Besides, one can see the presence of separate pikes for stars with 
A/Z = 3; 4; 5... and with A/Z = 3/2 on FigH 

It is important to underline that the measured mass of the Sun is in a 
good agreement with the claim obtained above and stating that the Sun core 
must be basically composed by hellium-10. This two measurement - the mass 
measurement and the frequencies measurement - build a test practically for 
the whole GIEP-theory. The first one tests the formula of mass and the other 
checks formulas for the radius of the core and the sound velocity (i.e. the core 
density). The agreement of this measurement results confirms the reliability of 
the obtained formulas and of the whole approach. 
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6.2 The star magnetic fields 



Another effect which follows from GIEP existence is the generation of a magnetic 
field by celestial bodies. These bodies, as they have electrically polarized cores, 
must induce magnetic moments due to their rotation: 

fita , (27) 

c 

where fi is the rotation velocity, Qbond = ^fPbondR 3 - 
Since the angular momentum of a star 

L re MQR 2 , (28) 

we can conclude that the giromagnetic ratios of celestial bodies must be directly 
expressed through world constants 

tf=^~^ (29) 
L c 

This theoretical predication can be checked by comparison with the mea- 
surement data. The values of giro-magnetic ratios for all celestial bodies (for 
which they are known today) are shown in FigEl The data for planets are 
taken from [5], the data for stars are taken from ^U], and those for pulsars - 
from Therefore, for all celestial bodies - planets and their satellites, Ap- 
stars and several pulsars, which angular momenta distinguish within more than 
20 orders - the calculated values of the gyromagnetic ratio fEo.l|29jl) agree with 
measurements quite satisfactorily with the logarithmic accuracy. 



7 Conclusion 

First conclusion, which we obtain from the above analysis states that there are 
four distributions, obtained from measurements, which depend on properties of 
the substance inside stars and which must be explained theoretically. The astro- 
physical models which are based on the Euler equation can neither explain the 
dependence of the velocity of periastron rotation from a chemical composition 
(Fig.0J), nor the star mass distribution (Fig.Q), nor their giromagnetic ratios 
(Fig.©). In this way, the quantitative agreement can be obtained only by an 
individual fitting of model parameters. The explanation of the spectrum of the 
solar oscillations (Fig.®) by means of series expansion on spherical harmonics 
(Ea.l|16|l) can be considered as a fitting only because its parameters Ai/q, £o 
and Dq are free and can not be obtained from the theory. 

Quite the contrary, taking into account GIEP effect opens possibilities for the 
quantitative explanation with acceptable accuracy of all measured data without 
using any fitting parameters. A good agreement of the relatively simple formulas 
and the measurement data has the easy explanation: the cores of stars consisting 
of hot dense plasma are well described by the known ideal gas formulas with 
small corrections, which are also well determined by modern plasma physics. 
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Figure 5: The observed values of the magnetic moments of celestial bodies vs. 
their angular momenta. On the ordinate, the logarithm of the magnetic moment 
over Gs ■ cm 3 is plotted; on the abscissa the logarithm of the angular momentum 
over erg ■ s is shown. The solid line illustrates Ea. H29(l . The dash-dotted line is 
the fitting of the observed values. 
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